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SCATTERING  THEORY  FOR  THE  ACOLYTIC  WAVE  EOUATION 
IN  AN  ARBITRARY  EXTERIOR  DOMAIN 


INTRODUCTION 

The  following  is  a preliminary  report  on  some  recent  theoretical  investigations  per- 
taining to  scattering  theory.  Further  study  and  the  physical  implications  of  these  results 
will  be  discussed  elsewhere.  The  scheme  for  direct  scattering  theory  is  outlined  first,  and 
an  important  topic  is  indicated  for  the  inverse  problem. 

If  u(x)  is  the  difference  between  the  instantaneous  pressure  and  the  equilibrium 
pressure,  p(.r)  is  the  equilibrium  density  of  the  medium,  c(x)  is  the  local  speed  of  sound, 
D.  is  an  open  connected  subset  of  Rn  with  bounded  complement,  D,  denotes  3/3/,  and 
a = c2(jc)p(.r)V  *l/p(v)V,  then  the  acoustic  wave  equation  is 

Dfu  = - Au,  t€R,  jell,  (1) 


wish  the  initial  conditions 

u(d,.t)  = fix)  and  DfU(O^c)  =g{xt,  (2) 

and  the  generalized  Neumann  bo<:ndarx  condition 

f c'2(r)p'Vx)Au(.vRlv|  dx  » f p"‘  <.r ) u[x)  • 7v(x)  dx,  x CE  12  and  i\7e€  L'iil).  (3) 

If  the  domain  S2  considered  here  has  a smooth  boundary  3 H,  say  C2,  then  boundary  con- 
dition (3)  is  •'quivalent  to  the  classical  Neumann  boundary  condition 

v • .wi /p.)  = 0.  x <3  312.  (4) 

with  v denoting  the  outward  unit  normal  of  3S2  at  .v 

This  acoustic  wave  propagation  problem  is  considered  here  as  a perturbed  system  in 
contrast  to  the  following,  named  “the  unperturbed  system'  : 

Djvltjc\  = ■ -40c(/^v) 

= - 7 2r(f  jrl.  .v  c SI,  (5) 


Nolo  Manuscript  sub.'.uttrd  June  1J,  1076 


1 


WfMpiWK*'-!***!*.';! 


CHEN  AND  YANG 

with  the  same  initial  and  boundary  conditions  (2)  and  (J)  of  the  perturbed  problem. 

We  impose  here  the  general  assumptions  that  apply  throughout  thL>  report: 

Assumr  .ten  1.  The  exterior  domain  12  has  the  finite  tiling  properry:  there  exist  an  open 
set  G in  Rn,  compact  sets  K 1 , ....  KN  in  Rn,  and  nonzero  vectors  xd),  . . . , jcW  such  that 


and 


9f2  C 0, 

0012  C UKj,  Kj<N, 

# 

{x  - *0  + tx^-.  o < / < i } c 12,  p.0  e nr\Kj. 


(6) 

(7) 

(8) 


This  pioperty  of  an  exterior  domain  is  due  to  Wilcox  (1].  Here  we  would  not  ex- 
clude the  case  12  ~ Rn. 

Assumption  2.  The  density  function  p(x\  is  C2(12)  and  real  valued,  and  for  some  constant 
J>  1 


j > pu,  > , x 12. 


Also, 


and 


p(x)  — 1 when  1 x i -*•  00 


(9) 


(10) 


r behaves  like  o{  lx  T1 ) when  Ixl-*'®® 

PU) 


(ID 


for  all  a,  1 < iai  = + . . . + an  *5  2. 

Assumption  3.  The  local  speed  c(x ) is  C'(12)  and  real  valued,  and  for  .nine  constant 

K > 1 


Also, 


rnd 


K & c'x ) > K"1 , x e 12. 


c(x)  -k  1 when  Ixl  -'  «> 


V In  c(x)  = o(!xTl . when  Ixi  -*  «». 


(12) 


(13) 


(14) 


Assumption  4.  The  “Stummel  condition”  is  satisfied  by  q(x)  = c'2tx)p*’ (x),  that  is,  for 
some  a > 0, 
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sup  f lq(y)  I2  lx-y  rn+4_a  dy  < + if  n > 4,  (15) 

where  the  integration  variable  y runs  in  the  disk  [ i*-y  I < 1 } nft  and  the  supremum  is 
taken  on  x £ Q. 

The  differential  operator  AQ  defined  by  (5)  on  ft,  with  assumption  1,  subject  to  the 
generalized  Neumann  boundary  condition,  is  a self  adjoint,  nonnegative  operator  in  the 
Hilbert  space  L2(ft).  Its  spectrum  is  the  closed  interval  (0,  °°)  and  is  absolutely  (spectral) 
continuous  and  without  eigenvalues.  These  interesting  results  are  proved  by  Wilcox  (1). 
With  these  four  general  assumptions  we  will  show  that  the  operator  A defined  by  (1) 
subject  to  the  generalized  Neumann  boundary  condition  is  also  a self-adjoint,  nonnegative 
operator  in  the  same  H'lbert  space  L 2(ft)  and  that  its  spectrum  contains  the  interval 
(0,  °°),  is  contained  in  (0,  <*>),  and  is  absolutely  continuous.  The  only  uncertainty  occurs 
when  the  origin  is  to  be  an  eigenvalue.  These  results  are  discussed  in  the  next  section. 

The  scattering  operator  S is  unitary  if  the  Miller  wave  operators  Wt  are  orthogonal, 

W'Wt  = /,  (16) 

ar.d  are  complete, 

W.  W*  =/  -£(0+),  (17) 

where  EQ.)  is  the  resolution  of  the  identity  for  A.  Therefore,  in  the  third  sect’on,  we  will 
discuss  the  existence  of  the  Miller  wave  operators  W.  and  properties  (16)  and  (17). 

7he  inverse  problem  of  the  scattering  theory  is  to  construct  c(x)  and  p(x)  from  the 
“scattering  amplitude.”  The  study  of  the  one-dimensional  piasma  inverse  problem  by  Szu, 
Yang,  Ahn,  and  Carroll  j:?l  and  the  articles  by  Faddeev  (3)  and  Newton  [4]  on  .he 
three-dimensional  inverse  scattering  problem  for  tie  Schrodinger  equation  all  center  on 
discussion  of  the  integra’  equation 

u{xjr,w)  = eix’k  +j'G(x-yrkjv)A'u(y,kjv)  dy  (iS» 

for  all  y such  that  yw  > x*w  for  each  unit  vector  w in  R,  where  A'  - A - A-  and 
where  G(x,fc;u’)  is  a certain  class  of  Green  functions  of  the  reduced  equation  from  the 
unperturbed  system: 


(A0  + k2)G(x,k)  =6(x).  (19) 

This  particular  problem  is  connected  with  the  direct  problem  and  is  discussed  in  the 
fourth  section. 

Other  references  relating  to  the  results  can  be  obtained  through  those  cited  in  this 
report. 
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SOLUTIONS  OF  THE  ACOUSTIC  WAVE  EQUATION  WITH  THE  NEUMANN 
BOUNDARY  CONDITION  AND  ITS  SPECTRAL  THEORY 

A solution  of  the  mixed  initial-boundary-value  problem  of  Eqs.  (1)  through  (3)  will 
be  constructed  by  using  a famous  spectral  theorem.  Then  the  absolute  continuity  of  the 
spectrum  of  A will  be  is  studied.  The  resulting  theory  provides  a preparation  for 
constructing  the  wave  ot  orator  and  scattering  operator. 

We  recall  here  the  initial-boundary-value  problem: 

Dfu  = - Au,  t > 0, 

= c2(x)p(x)7‘Hp\x)Vu,  x €12,  (20) 

u(0,  x)  - fix)  and  £)tu(0,  x)  =g(x),  xGfi;  (21) 

fc'2(x)p~*(x) Au(:<)v(x)  dx  = | p~l(x)7u(x)'Vv(x)  dx,  x G 12  and  i;,7o  G L2(Sl).  (22) 


The  formulation  of  the  problem  will  based  on  the  following  function  spaces: 


u G L2(f2)  **j  1 u(x)  i2  dx  < <*>,  x G £2,  in  the  Lebesque  Pleasure, 

(23) 

u e **  Dau€L2{Sl),  lor  1 < m. 

(24) 

u G L2(ft;c2p)  ( c2p)~xl2u  G 

(25) 

u G Hm(Sl;c2p)  (c2."!'1 ,2rnu  G L2(ft),  1 ^ m> 

(26) 

and 

u G Xj  **  p*l2VJu  G L2(12). 

(27) 

These  spaces  are  Hilbert  spaces  with  respect,  to  the  following  inner  products  respectively: 

(u.i  ) = I*  u(x)vix)  dx,  x € SI, 

(28) 

= £(0°u,  D°o),  iaKm, 

(29) 

< u,v ) = fc'2(x)p*J(x)  u(x)  v(x)  dx,  x G fl, 

(30) 

(u.i>)r.;  = £ <£°u,  D°f>,  >al  < m. 

(31) 

[u,ul,  = y f VT:(x vfx)  dx,  x G ft  and  1 < / < i < 2. 

p(x) 

(32) 
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It  is  clear  from  the  characteristics  (9)  and  (12)  of  p(x)  and  c(x ) that  L2(£l)  - H0(11) 
(which  is  equivalent  to  L2( i?-;  c2p)  - //0(f2;  c2p),  in  the  view  of  their  no.-ms). 

Suppose  that  the  boundary  dS2  of  n is  sufficiently  smooth,  say  C2.  Then  Green’s 
theorem  implies  that 


■k 


Vu(x) 


v(x)  dx 


'P(x) 
■/['  • 


7<(x)  • Vvfx)  dx 


— • jpu(x)  u(x)  ds(x),  x e an, 
P(X  J 


(33) 


where  v is  the  unit  outward  normal  of  d!2  at  x.  This  means  that 

-<Au,v>  + {u.t/J,  = /*■"■[»'  • 7u(x)J  * u(x)dS(x),  xEdtt. 
J P(x) 


(34) 


Therefore,  for  u € X2,  u satisfies  the  classical  Neumann  boundary  condition  v * 7 u(t,x) 
0 if  and  only  if  v satisfies  the  relation  (22),  which  in  the  new  notation  is 


<Au,w>  = v € XxCiL2(Sl\c2p) 


(35) 


Definition  1.  A function  u € X2  is  said  to  satisfy  the  generalized  Neumann  condition  if 
and  only  if  (35)  holds. 

This  definition  does  not  require  the  assumption  of  the  smoothness  of  the  boundary 
dn  of  n.  However  it  defines  the  classical  Neumann  condition  if  dfi  is  smooth. 

Furthermore  a definition  is  introduced  for  a closed  subspaco  m X2: 

u € 50  «*  u G X2  and  satisfies  (35).  (36) 

The  construction  of  a solution  of  the  initial-boundary-value  problem  is  based  on  the 
linear  operator  A in  H2  given  by  (20)  v. : domain 

D{A)  = X.  (37) 

Theorem  l.  A is  a self-adjoint  operator  on  the  Hilberi  space  L2(Slx2p).  Moreover  A ? 0. 

The  verification  of  the  assertions  is  based  on  the  following  result: 

Lr  nma  1.  Let  X be  a Hilbert  space  and  let  L:  X -+X  be  a linear  operator  dentely  de- 
fined in  X.  Assume  that  L C L*,  the  adjoint  of  L,  L > 0.  and  that  the  range  R(l  + L) 
of  l + L is  X.  Tht-n  L is  self-c/ijc-int. 

Proof  of  Lemn -.  1.  L > 0 indicates  that  the  deficiency  indices  of  the  symmetric  operat  >r 
L and  of  £»*,  an  extension  of  L,  ere  equal  (5,p.268|.  Condition  R(l+L)  - X implies  the 
deficiency  index  is  zero  i-'k*.  hence  L is  self-adjoint. 
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Proof  of  Theorem  1.  Verify  the  conditions  of  L>mma  1 for  K * L2(H;c2p)  and  L = A. 
The  space  CQ(Sl)  of  infinitely-many-times  continuously  differentiable  functions  with  com- 
pact support  in  $2  is  a subset  of  D(A)  = Jf;  hence  D(A)  is  dense  in  Let  u and 

v be  any  two  elements  of  D(A).  Then  (35)  and  (32)  yield 


C4u,  v)  * [u,  u]  i = f~- - ?u(x)  ♦ Vt>(x)  dx  = (u,Av).  (38) 

J r\*l 

Consequently  D(A)  C D(A*)  and  A*u  = Au  for  all  u € D(A);  that  is,  A C A*. 

The  assertion  A > 0,  also  one  of  the  conditions  to  be  checked,  follows  from  the 
result  yielded  by  (35)  and  (32)  that  for  all  u € D(A) 

<Au,u)  = l«.«h  =j\p'll2(x)Vu(x)\2dx,  xen.  (39) 

The  only  condition  left  to  be  verified,  R(I+A)  = L2(fi;c2p),  means  that  for  each  f 
in  L2(fl;c2p)  there  exists  an  element  u of  D{A)  such  that 

(u,v)  + (Au , v ) - (f,  v > for  all  v in  H1(n#2p).  (40) 

This,  together  with  (35),  is  equivalent  t'* 

( u , v > + [u,  uj  j = <f,  v)  for  all  v in  //t(i2;c£p).  (41) 

Sufficiently,  if  the  equivalent  inner  product  tor  Hx(Sljc2p)  that  we  use  is 

<u,  v)  + u,v  G Hx(Sl;c2p),  (42) 

we  need  to  verify  the  existence  of  u such  that 

{u,  v } = (f,  v\  v€Hx  (fl;c2p).  (43) 

However, 


|</.  v)l  < (f,ft*l2(v,v)112  < const  {f,f}l,2{ v,  v}112. 

The  Riesz  representation  theorem  in  the  Hilbert  spac.  (Hx (fl;c2p),  {. , .})  yields  the 
existence  of  an  element  u in  Hx{Sl;c2p)  satisfying  (43)  and  then  (41).  On  the  other  hand 
(41)  implies  (40)  for  all  v in  cj(n).  Thus  Au  * f - u with  the  member  on  the  right  side 
in  L2(I2,c2p);  hence  Au  is  also  in  L2(fi,c2p).  Moreover,  because  the  validity  of  (40)  itself, 
is  implied  when,  Cq(S1)  is  dense  in  f/j(fl,c2p),  the  combination  of  (40)  and  (41)  ensures 
that  u satisfies  the  generalized  Neumann  condition  '.35).  Thus  R(!+A)  = X is  verified,  and 
the  proof  of  theorem  1 is  complete. 

Therefore  the  Kato  |5,p.331)  second  representation  theorem  ensures  the  following 
corollary. 
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s Corollary  1.  A has  a nonnegative  square  root  A112  whose  domain  D(A 1/2 ) = //j(ft;cp1/2) 

\ has  the  inner  product 

l {u.v}  = 22  J c'x  (x)p~ll2(x)Dau(x)Dav(x)  dx, 

| Furthermore  A112  satisfies  the  relation 

\ 

? {Ail2u,All2u)  = 2 2(P~ll2(x)DjU,p~ll2DjU), 

\ 

/;  From  the  results  in  theorem  1 and  corollary  1,  the  argument  of  Wilcox  { 1 ) , with 

j slight  adjustment,  gives  the  following  theorem. 

♦ 

j Theorem  2.  For  each  f i.i  D(A)  and  g in  D(A 1/2)  there  exists  a uniquely  defined  strict 

l solution  u with  finite  energy  of  the  initial-boundary-value  problem  (20),  (21),  and  (35) 

with  t £ R such  that 

\ u e c2[b,  L2(fl,c2p) ) ncl[fl,tfi(*vp1/2))  nc(«,->;), 

l 

e and  u has  the  energy  integral  in  the  two  equivalent  forms 

5 

» 

l E(u,n,t)  = {Dtu(t),D,«(t)}  + 22  {Djn(t),Dju(t)),  j = 1,. . . ,n, 

« {Dtu(t),D, «(()}+  {Axl2u(t),A'l2u(t)\  (46) 

and  has  the  constancy  of  energy 

E(u,  ft,  t)  = Z iD,f-  *>jf}  * { g . g}.  (47) 

From  the  spectral  theorem  for  A and  the  associated  operator  calculus,  we  have  the 
following  theorem. 

Theorem  3.  For  real-valued  functions  f in  L2(ft)  and  g in  D(A'112),  define 

h = f + iA'l,2g  e f,2(ft).  (48) 

Then  the  solution  in  L2{$1)  defined  by 

u(t)  = (cos  tA*l2)f  + M*1/2  sin  M1/2)g  (49) 

with  bounded  coefficient  operators  satisfies 

U(t,x)  = Rev(t.x),  (50) 

where  e( t,  x ) is  the  complex-valued  solution  in  L2(ft)  defined  by 

t/»,0)  = e'ilA'll2h.  (51) 


x€ft,la|<l.  (44) 


l<;<n.  (45) 
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Because  the  operator  A is  nonnegative,  its  spectrum  o(A ) is  contained  in  the  interval 
(0,-  ).  The  nonexistence  of  the  positive  eigenvalue  is  ensured  by  the  following  theorem. 

Theorem  4 (Mochizuki  [6].  Assume  assumptions  2 through  4.  Then,  u - 0 in  ft  is  the 
only  L2  {&)  solution  of  the  equation 

V • Vu  + Xc'2<x)p‘1(x)u  = 0,  x€ft.  (52) 

Proof.  It  suffices  to  check  the  conditions  imposed  by  Mochizuki.  His  first  three  condi- 
tions are  presented  by  assumptions  2 through  4 by  setting  a;fc(x)  = y/jk  p_1(x),  6-  = 

0,  and  q(x)  = Xc*2(x)p‘1(x).  Particularly,  if  n < 3,  the  “Stummel  condition” 

supj  lq(y)|2  dy  < °°,  (y  € ft:  lx-y|  < 1,  x 6 ft)  (53) 

is  a consequence  of  the  boundedness  and  the  smoothness  of  c(.r)  and  p(x)  on  ft.  His 
fourth  condition  is  given  by  his  remark  1.2  and  by  (9),  (11),  and  (14).  As  for  the  unique 
continuation  property  in  his  last  condition,  the  smoothness  assumption  on  c,x)  and  p(x) 
implies  the  Holder  condition  ( 7,  p.  23],  which  gives  the  property  (according  to  his  remark 
1.1,  or  also  according  to  Ref.  8 or  Ref.  9). 

Subject  to  some  adjustment,  the  employment  of  the  standard  principle  of  limiting 
absorption  yields  the  absolute  continuity  of  the  continuous  spectrum  of  A.  The  detailed 
proof  will  appear  in  a follow-up  article.  The  precise  statement  is  the  following. 

Theorem  5.  The  resolution  £(*)  of  the  identity  for  the  operator  A is  absolutely  continu- 
ous on  any  closed  interval  in  (0,  «>). 

Remark  1.  Whether  X = 0 is  an  eigenvalue  or  belong  to  the  continuous  spectrum  is  not 
yet  clear. 


DISCUSSION  OF  THE  M0LLER  WAVE  OPERATORS 

As  indicated  in  the  Introduction,  all  properties  of  A studied  in  the  last  section  hold 
for  Aq.  Moreover  the  spectrum  of  A0  is  {0,°°)  and  is  absolutely  continuous.  These  are 
results  of  Wilcox  { 1 ] . These  results  and  those  of  the  last  section  ensure  the  existence  of 
the  unitary  groups  and  e~iLt,  < t < °°,  associated  with  the  self-adjoint  operators 

Lq  = A}/2  and  L - A'12. 

The  strong  limits 

Wt  = S - lim  e,tAe',,'°t . t -*  (54) 

are  called  the  Mdller  wave  operators.  This  yields 

lim|le',*'</i-e',i'0,htii  = 0,  t -*  ♦ <».  (55) 


and 
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h = Wt  h±.  (56) 

The  following  map  is  the  scattering  operator: 

(57) 

It  is  required  to  be  unitary  on  L2(£l),  which  is  an  ea..,  consequence  of  the  orthogonality 
(16)  and  the  completeness  (17)  of  the  Miller  wave  operators. 

The  existence  of  the  strong  limit  (54)  is  proved  by  employing  the  vanishing,  when 
t -*■  t °o,  of  the  local  energy  of  the  solution  to  the  unperturbed  system  and  employing 
the  decaying  in  rate  -1  - 5 of  lt|  for  the  first  second-order  derivative  of  the  solution  to 
the  unperturbed  system. 

The  coincidence  of  the  Miller  wave  open  tor  W+  with  the  stationary  wave  operators 
U±  and  the  properties  of  the  orthogonality  and  the  completeness  for  t/±  guarantee  the  cor- 
responding properties  for  W± . The  proof  of  the  coincidence  and  of  these  two  properties  for 
Ut  is  based  on  the  expansion  principle  of  the  generalized  eigenfunctions  of  A,  which  has 
been  employed  by  a dozen  different  authors,  including  Ikebe  [10] , Mochizuki  [11,12] , and 
Wilcox  [1] . The  details  will  be  given  in  the  follow-up  article.  A portion  of  the  article 
involve*  the  principle  of  limiting  absorption,  which  is  studied  in  the  next  section. 


PRINCIPLE  OF  LIMITING  ABSORPTION  AND  THE  INVERSE  PROBLEM 

There  is  no  doubt  about  the  role  played  by  the  principle  of  limiting  absorption  in 
the  eigenfunction  expansion,  which  is  a keystone  for  scattering  theory.  However,  in  the 
three-dimensional  inverse  scattering  theory  for  the  Schrodinger  equation,  both  Faddeev 
[3]  and  Newton  [4]  studied  the  “Volte rra”  integral  equation,  where  the  principle  of 
limiting  absorption  again  played  an  essential  role.  We  will  apply  the  principle  of  limiting 
absorption  to  the  acoustic  wave  equation.  This  should  bring  some  light  to  the  study  of 
the  inverse  problem  of  the  scattering  theory  for  the  acoustic  wave  in  three  dimensions, 
which  will  extend  Ref.  2. 

In  the  investigation  cr  recovering  the  perturbed  operator  A1  = A - A0,  the  reduced 
wave  equation  of  acoustics  without  an  obstacle, 

Au[x,k)  + k2u(x,k)  = 0,  x, kE.R3,  (58) 

is  expected  to  have  a set  of  solutions  u{x,k;w)  such  that  the  function  u(x,  k;w)e~'k’x 
has  an  analytic  continuation  into  an  upper  halfplane  in  the  variable  s = ft1  • 7 at  a fixed  x 
and  kr  = k - (fe*7)7;  and  u(x,k;w)e~i"’x  - 1 decreases  with  large  |s|.  From  an  alternative 
form  of  (58),  namely, 

AqU  + k2u  = A jU 


= fc2[l  - c2(x)p2(x)]u  - 


VP(*) 

P(x) 


Vu. 


fe.x  € R3, 


(59) 
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we  should  study  the  integral  equation 

u(x,  ft;  w)  = e'k'x  + I G(x  - y,  ft;  w)A , u(y,  ft;  w)  dy  (60) 

with  x,  ft  in  R3  and  with  a unit  vector  w,  where  the  Green  function  G(x , ft;  u;)  is  a 
limiting  value  of  the  following  function  with  Im  s = 0: 

G{x , y;  s)  = (?n)'3J  e‘^m*sw)‘x[(m  + sw )2  - y2  - s2)"1  dm,  (61) 

with  m G R3,  w = ql\q\,  s = k‘w  + t'h/l,  and  y2  = ft  - ( k‘w )2.  As  indicated  by  Faddeev, 
the  irregularities  in  the  expression  under  the  integral  sign  are  proven  by  the  equations 

m2  = y2  and  m-w  = 0 (62) 

and  do  not  depend  on  s with  Im  s t 0.  This  implies,  through  a singular  Fourier  integral, 
that  the  function  defined  by  (61)  is  an  analytical  function  of  s in  the  upper  half  plane 
and  that 


|G(a,  y,  s)  exp  [Im  s(x-u>)]  i < const  f(l  + Ixl)*1.  (63) 

The  existence  of  an  analytical  continuation  of  the  Green  function  G(x,  ft;  w ) in  the  upper 
halfplane  in  the  variable  s = k'W,  which  satisfies  (63),  is  just  the  characteristic  property 
that  defines  this  function  uniquely. 

On  the  other  hand,  in  (59)  is  still  a differential  operator.  We  remark  that 

Vpl  y) 

G(x  - y.  ft;  w)Vu(y , ft ; u;)- 

p(y) 


* G(x  - y,  ft;  ie) 


Vu(f,  k;  w) 


Vp(y) 

P(x) 


- u(y,  ft ; w) 


Vp(x) 

. P(*) 


\Vp(x)\ 

P(x) 


21 


= V, 


G(x  - v,  ft;  u?)u(y,  ft;  w) 


?P(y)' 

P(.v) 


\’x.yG{x  - y,  ft;  u>) 


yp(y) 

p(y) 


?2p(s) 
p(x ) 


|?P(j:)12 

P(.v) 


G(.r  - v;  ft,  w ) 


u(y,  ft;  u>). 


(64) 


Formally  this  relation  and  the  divergence  theorem  of  Gauss  give  another  representation 
of  (60): 

u(: x,  ft;  it’)  = e‘h'x  + ) G(x  - y,  ft;  u>).42u(y,  ft;  re)  <fv 

+ f'.x..  G(xpy,  ft:  «•)•  ' P(-  > u(v,  ft;  ir)  dy.  (65) 

J T ' p(y) 
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where  the  multiplication  operator  A 2 is  defined  by 


A2u  = d2[  1 - c2(y)p2(y)} 


'V2P(  y) 
. P(y) 


lvp(y)!2' 
p(y)  . 


u(y,  fe;  w). 


Here  the  following  condition  is  assumed: 


lim  J G(w 


VP(v) 

y,  k;  u>)u(y,  k\  w)y  • dS(y ) = 0,  iyl  = r -*•  », 

p(y) 


{66) 


(67) 


with  the  surface  integral  over  the  sphere  with  radius  r.  However  some  u(y,  k;  w)  is  re- 
quired to  be  bounded,  and  since  G(w  - y,  k ; u>)  has  the  asymptotic  estimate  (63),  condi- 
tion (11)  in  assumption  2 ensures  (67).  The  representation  (65)  is  not  justified  yet, 
because  the  asymptotic  estimate  for  VyG(y,  k ; w ) at  large  |y|  is  not  clear. 

This  brings  in  the  issue  concerning  the  asymptotic  estimate  of  both  G(y,  k;  w ) and 
VyG(y,  k;  w)  given  by  (61)  with  the  principle  of  limiting  absorption  employed. 

In  the  rest  of  this  section,  we  will  follow  the  idea  used  in  Ref.  13.  First  we  recall  a 
well-known  formula  [11,  13] : 


£ 


e,rp[r  - (a  + »6)]‘ 


dr  = nr(sign  b)H(p  sign  b)ei^a*,b)p. 


(68) 


where  H(t)  is  the  Heaviside  function.  From  (61)  we  have  for  a multiple  index  a 

D«G(x , p;  s)  = (2nr3eigw'xJ  r2<7+  + J.)  dr. 


(69) 


and 


7*  + J.  = j ^ [i(mr  + $ic)lu(p2  «•  (sm  •w)2]~1,2e'rm  + (snrw) 

+[p2  + (sin-  u/)z]  1/2|  + |r  + (sin*  u;)  - [p2 

+ (*m- w)2 )1/2}  jdS(m). 


(70) 


Without  loss  of  generality  we  can  assume  ffrx  - x^m j.  By  the  Morse  transformation  [14] 
we  choose  a local  coordinate  (r;,  ,7j2),at  (±  1,0,0)  si  "h  that  m x = ± (1  - t;2  - q2).  Then 
the  estimate  for  large  |x|  is 
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«/+  + J.  = const  ~ ( p2  + (sx*u>)2]'1/2(e,r**l[,'(*r  + su;)]Q({r  + (sx'w) 

+ [p2  + (sx'u;)2]^2}'1  + {r  + (sx*u;)  - [p2  + (ac -w)2  ] 1/2)_1 } 

+ c',r'*l[i(-  xr  + su/)]a({r  - (sx'w)  + [p2  + (s5c  *u?)2  ] 1/2 }_1 

+ {r  - ( sx'W ) - [n2  + (sx*u/)2]1/2}'1)]  + q(rx),  (71) 

where  q(x)  is  estimated  for  large  |x|  as 

lq(*)i  + l\/g(x)|  < const  |x|*2.  (72) 

Because  the  denominators  in  (70)  would  not  vanish  when  m*u;40  (by  (62)),  there  is  no 
singularity  in  (71)  when  5-ur+O.  After  substitution  of  (71)  and  (72)  into  (69)  and  then 
integration  by  parts,  we  have  the  estimation 

D^G(x,  p;  s)  < const  |x|'2,  if  x'wtO.  (73) 

Assume  now  that  x'w  = 0.  Then  (71)  becomes 

+ J.  = const  2ff(r|x|p)'1  { tfrix\fxrm)a  + e'ir^[i(sw  - xr)) a}[(r  + m)"1 

+ (r-m)_1]  + q(rx),  (if  3c -to  = 0).  (74) 

Let  0(r)  be  a CT  function  with  compact  support  such  that  0(r)  has  values  between  0 and 
1 and  is  equal  to  1 at  r = ± p.  Therefore  (69)  and  (74)  imply  that 

D“G(x,  p;  s)  = (27 r)~3j'  r2(J+  + «/_)[l  - p(r))  dr  (=  /t) 

+ (23r)‘2(lx|p)‘1  const j (e,r|x,p4.(r;  3c,  s,  w,  a) 

+ e*,r|jt|0.(r;  x,  s,  w,  <\M[r  + p)*1  + (r  - p)'1 ) dr(=  l2) 

+ 0(|xl‘2),  ifx-u>  = 0.  (73) 

where  <pt  - r^(r)[«(sa±xr)]a.  There  is  no  tingularity  in  the  integral  for  / j.  Integration  by 
parts  yields  / j = 0(ix|'2).  Employment  of  .68)  for  /2  gives 

l2  = (2ir)x]p)'1  const  [^+(|xl;  x,  s,  w,  o) 

+ v.(  - ixl;  x,  s,  w,  a))  * (sin  (mix  !)//( lx|)l,  (76) 


It 
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where  f * g means  the  convolution  of  f and  g and  where  the  function  y±(l*l;  x , s,  w,  a) 
is  the  Fourier  transform  of  <t>+(rjc,  s,  w,  a)  in  r.  Since  0+  is  compactly  supported  in  r, 

0+  is  a smooth  rapidly  decaying  function  in  |x|.  Therefore 

D%G(x,  n;  s)  = OKIxl/i)*'  ] , for  large  |x|  with  5 ‘w  = 0),  (77) 

for  any  order  a of  derivatives. 

Then  for  either  x-w  = 0 or  x-w  1 0,  we  have 


I D*G(x,  p:  s)  exp  {/m  s(x*u;)}l  < const  (1  + |/xxj)_1 , 0 |u|  < (78) 


which  is  (63)  when  a = 0. 

This  estimate  ensures  the  two  integrals  on  the  right  side  of  (65)  are  well-defined 
under  the  general  assumptions  for  c(x)  and  p(x).  Equation  (65)  will  be  further  studied  in 
the  follow-up  article. 
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